Flow profiling of a surface acoustic wave nanopump 



^1" 
O 

o 

(N 



Z. Guttenberg 1 , A. Rathgeber 1 , S. Keller 2 , J.O. Radler 2 , 
A. Wixforth 3 , M. Kostur 3 , M. Schindler 3 , P. Talkner 3 
1 Advalytix AG, Eugen-Sanger-Ring 53.0, D-85649 Brunnthal, Germany. 
2 Ludwig-Maximilians-Universitat, Geschwister-Scholl-Platz 1, D-80539 Miinchen, Germany. 
3 Institut fur Physik, Universitat Augsburg, Universitatsstrasse 1, D-86135 Augsburg. 

(Dated: February 2, 2008) 

The flow profile in a capillary gap and the pumping efficiency of an acoustic micropump employing 
Surface Acoustic Waves is investigated both experimentally and theoretically. Such ultrasonic sur- 
face waves on a piezoelectric substrate strongly couple to a thin liquid layer and generate an internal 
streaming within the fluid. Such acoustic streaming can be used for controlled agitation during, e.g., 
microarray hybridization. We use fluorescence correlation spectroscopy and fluorescence microscopy 
as complementary tools to investigate the resulting flow profile. The velocity was found to depend 
on the applied power somewhat weaker than linearly and to decrease fast with the distance from 
the ultrasound generator on the chip. 
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I. INTRODUCTION 

For a biological hybridization experiment, typically, 
single stranded DNA is covalently bound to small probe 
spots on the surface of a glass slide, forming a so-called 
micro array. A solution containing fluorescently labelled 
sample molecules is then spread across this micro array 
where the molecules are allowed to bind to the respective 
probe spots. Successful hybridization is usually observed 
as an increase in the fluorescence intensity on a spot. To 
have a reference, in co-hybridization experiments a two 
color coding is used. Usually, such experiments are per- 
formed in sealed chambers to prevent the thin liquid film 
above the slide from evaporation. Given the typical film 
thickness of less than 100 microns, the system exhibits 
a very small Reynolds number, and the molecules in the 
fluid can only move diffusively [lj. Diffusion, however, 
is a very slow process. Even for small molecules such 
as Rhodamine the distance travelled in water after 10 3 s 
is only about 1mm. For relatively large DNA molecules 
the diffusion driven displacement will be even smaller. 
On the contrary, the binding reaction of complementary 
DNA strands is a relatively fast process. Especially in 
the case when the concentration of probe molecules on 
a spot is much larger than the one of their counterparts 
in solution, most of the molecules close to their binding 
partners at the glass slide will become immobilized af- 
ter a short time. This can lead to the formation of a 
depletion zone with a lower concentration close to the 
solid surface. To reach the equilibrium concentration of 
bound DNA, molecules have to diffuse through the deple- 
tion zone, which becomes the time limiting bottleneck for 
the hybridization process. Typical hybridization assays 
are therefore performed over night or even over week- 
end, and then are stopped at some more or less arbi- 
trary point. Such experiments are also referred to as 
"end point" experiments. Apart from the above men- 
tioned non-equilibrium conditions, shorter hybridisation 
times also lead to a lower fluorescence signal intensities. 
Quite obviously, this process can be accelerated if the 
molecules were to be actively transported along the solid 



surface plane. Due to the low diffusion constant of large 
DNA strands even a slow movement can speed up the hy- 
bridisation considerably. Another aspect is that hybridis- 
ation is often conducted with more than one kind of DNA 
with different fluorescent labels simultaneously. The flu- 
orescence intensities of the different immobilized sample 
molecules are then evaluated quantitatively to interpret 
the results. However, when the DNA binding reactions 
were not in equilibrium at the end of the hybridisation 
time, the results can be misleading. 

As the amount of DNA available for a hybridisation 
experiment is often very limited one has to reduce the 
volume as much as possible. Therefore, pump driven 
flow systems usually having large dead volumes are 
not convenient to induce the molecule transport during 
the hybridisation process. There are numerous other pos- 
sibilities to actuate thin fluid films, amongst them for 
example electrokinetic forces Q 0- 

In this report, we describe a novel technique to actuate 
a thin liquid film in a capillary gap. We employ surface 
acoustic waves (SAW) propagating on the surface of a 
piezoelectric solid to induce acoustic streaming and mix- 
ing within the fluid. To understand the fluid dynamics 
in such a system, we concentrate on how the acoustic 
streaming transforms into a velocity field profile within 
a homogeneous DNA solution in the capillary gap. 

A surface acoustic wave (SAW) can be excited by the 
application of a radio frequency signal (RF) to an in- 
terdigital transducer (IDT) on a piezoelectric substrate. 
The wavelength of the excited SAW is defined by the ge- 
ometry of the IDT. Typical displacement amplitudes of 
a SAW arc in the nanometer range, depending on the 
applied RF power. If the SAW supporting surface of the 
substrate is in contact with a liquid, and if the SAW 
has any displacement component normal to the surface, 
ultrasonic SAW power is leaking into the fluid in form 
of an acoustic wave. Therefore, this so-called leaky SAW 
(LSAW) is expected to decay exponentially with distance 
from its source. 

The interaction between the SAW and the liquid film 
leads to an internal, acoustically induced streaming. The 
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FIG. 1: The measurement setup for the evolution of dye spots 
consists of a 200 /im slit between two plexiglas plates (3x3 
cm 2 ) kept in the right distance with spacers. Integrated in the 
lower plate is the SAW chip, being connected to the RF source 
outside the fluid film. The transducer is located in the middle 
of the plate. Images were taken with a video microscope from 
above through the upper plate. 



FIG. 2: The graph shows the sketch of the 4x23 mm 2 LiNb0 3 
chip with gold structures. Situated on the right side is the 
interdigitated transducer(lxl mm 2 ) inducing the SAW per- 
pendicular to the long side of the chip. The contact wires are 
connected with the Rf power source on the left end. 




FIG. 3: Flow pattern generated by SAW. The fluid is pumped 
perpendicular to the metal fingers of the IDT and is sucked 
in parallel to them. 



reason for this is a net pressure gradient in the direction 
of the sound propagation in the fluid. Usually, this leads 
to the generation of a fluidic jet, which can be nicely 
visualized employing a dye solution and a fluorescence 
video microscope. The captured images are evaluated in 
a subsequent image processing step, being a versatile tool 
for the interpretation of dynamic processes. This method 
yields real images of the observed sample as a whole, and 
therefore reduces the problem of a false interpretation of 
the data. For more precise measurement of the flow field 
and velocity, we use fluorescence correlation spectroscopy 
(FCS) 5]. In this method, the movement of fluorescent 
particles can be determined in the focal spot of a confo- 
cal microscope by evaluating the autocorrelation function 
of the intensity fluctuations. This represents the rate of 
particles entering and leaving the small observation vol- 
ume. 



II. THE EXPERIMENTAL SETUP AND 
MICROFLUIDIC CHARACTERISTICS 

We assembled a special microfluidic arrangement to 
induce streaming in a narrow capillary gap between a glas 
slide and a plexiglass block, containing a SAW device. 
Fig. n gives a sketch of the principal setup as seen from 
the side. Fig. shows the chip architecture. The width 
of the gap is controlled by spacers, it can be continuously 
adjusted down to 30/im. In the measurements described, 
we worked with a spacers of 200 ^m. The SAW chip 
consists of a 128° rotated Y-cut, X propagating LiNbC>3 
single crystal chip. 

The IDT consists of 22 finger pairs, 0.85 mm aper- 
ture, and a period A = 28/im. The IDT fingers are ther- 



mally evaporated 300 nm thick Au electrodes, fabricated 
by a standard photolithografic technique. The complete 
chip was protected with a 800 nm thick, sputtered Silicon 
Dioxide layer, which was removed at the contact pads to 
allow for electrical connection. The IDT was driven with 
an AC signal at the resonance frequency /o, being de- 
fined as f — c/A. For a SAW velocity c rj 3800 m/s, 
the frequency fo is hence about 136 MHz. Such an IDT 
efficiently converts the applied RF signal into an acous- 
tic wave, which in this case is launched as a bidirectional 
beam perpendicular to the fingers. 

If the capillary gap is filled with a fluid, switching on 
the RF power leads to a flow pattern as sketched in Fig. [3] 
The fluid is expelled from the transducer region to the 
left and right side and flows back to the IDT from bottom 
and top. According to the high viscous damping the flow 
velocity decreases rapidly with distance from the IDT. 

Further investigations were focused on the character- 
ization of the flow field. The issue of this work was to 
find out how the setup can be optimized for agitating a 
large hybridisation area. For this purpose we first want to 
discuss the nature of the observed SAW-induced motion. 

The propagating LSAW couples to a sound wave in 
the liquid. Because the sound velocity in fluids cii qu id 
is always smaller than the LSAW-velocity in the solid 
substrate, sound is radiated into the liquid under an angle 
obeying Snellius' law of diffraction: 

• In \ ^liquid /-, \ 

sin(0fl) = —2 (1) 

clsaw 

This can be understood as a consequence of phase match- 
ing between the LSAW and the radiated sound beam in 
the liquid. This " mode conversion" leads to an exponen- 
tial decay of the LSAW. In our case the observed damp- 
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ing in power was about ldBm per wavelength, which is 
in good accordance with values found in refs. |(j and 
7]. This results in a decrease of 35dBm per mm. In the 
following we will give a short overview how the sound 
radiated in the fluid can induce a streaming force. 

Sound travelling through a liquid is attenuated by the 
viscosity along its transmission through the medium. If 
the intensity is high, this attenuation creates an acous- 
tic pressure gradient along the propagation of the wave. 
The gradient induces a force in the same direction that 
causes a flow in the fluid. This conversion of an attenu- 
ated sound wave into a steady flow is a nonlinear effect 
which is known as acoustic streaming. When a SAW, i.e. 
a surface acoustic wave with a spatially constant ampli- 
tude, is in contact with a fluid the acoustic streaming is 
only a minor effect 0. Much stronger is the streaming 
that results from the exponential decay of the LSAW. 
This decay causes the conversion into sound waves in the 
liquid within a short distance leading to a large gradi- 
ent of the sound amplitude and consequently to a strong 
force that generates the acoustic streaming. 

Here we do not want to go into the details of the 
rather complex mechanism of acoustic streaming @ j @ • 
We only mention that the effective force that causes the 
steady flow can be expressed in terms of the sound ve- 
locity field in the fluid 0: 



p((v f -V)v/+v/(V-v/)) 



(2) 



where the angular brackets denote time average that re- 
moves the fast oscilations of the sound wave. 

Finally we note that the force is proportional to lu 2 . 
High force densities therefore can be achieved by means 
of high LSAW frequencies. On the other hand, lower 
frequencies lead to smaller damping rates of the LSAW 
and the force density acts on a larger volume. For an 
optimal design a trade off has to be found. 



III. MATERIALS AND METHODS 
A. Fluorescence Microscopy 

a. Microscope To visualize the activity of the SAW 
transducer, a Na-fluorecein solution (lmg/ml) is pipetted 
into a small hole in the cover of the fluid film. For the ex- 
citation of the fluorescence of the dye, blue LEDs are at- 
tached to a ring some distance above the cover. Fluores- 
cence images are taken by a video microscope (Olympus, 
Hamburg) with a suitable filter to suppress the excitation 
light. 

b. Image processing Fluorescence images of the 
evolving Na-fluorescein spot (see Figs. 01 were taken 
every 10s and stored into an avi movie file. The avi movie 
was further processed by a procedure converting the pic- 
tures into binary images with the help of the Shannon 
entropy function and detecting the edges of the dye spot 
with the particle finder algorithm. The coordinates of ev- 
ery edge line were put together in one file (Fig.^J and the 



positions of the leading edge to the left and the right was 
extracted for every recorded point of time. This leads 
to a distance vs time diagram (Fig. 0) . Differentiating 
this data yields a speed vs distance diagram. Also, the 
area of the Na-fluorescein spot and the speed of the area 
increase can be calculated. 



B. Fluorescence Correlation Spectroscopy 

In 1972 Magde, Elson and Webb published the first 
work on FCS [l(j ■ In the following years the same authors 
elaborated the technique and it's applications in detail 
[ll| [l2| |l3j . where in the last reference theory and first 
experiments with FCS on systems with laminar flow are 
presented. Until the early 1990s FCS was mainly used to 
determine diffusion properties and kinetic constants of 
chemical reactions in solution. In 1993 the introduction 
of a confocal detection optics by Rigler et al. Q gave 
new impulses to the method. Schwille [l5|. Krichevsky 
|16| and Thompson |T^ have given reviews about the 
principles, applications and potentials of FCS. 

The raw signal in FCS experiments is the time- 
dependent intensity signal coming from a small and fixed 
volume (« 0.5 fl) created by a laser focus. Fluorescently 
labelled molecules entering and leaving the open volume 
lead to fluctuations in the recorded signal F(t): 



F{t) = qQ [ I(f)C(f,t)dV 
Jv 



(3) 



Here, q stands for constant intensity losses due to the 
experimental setup e.g. detector efficiency, filters. Q is 
the fluorescence quantum yield of the dye and C (r, t) de- 
notes the concentration of the labelled particles, 1(f) = 
Io exp(— 2(x 2 + y 2 )/ro) exp(— 2z 2 /z n ) is a combination of 
the collection efficiency and the excitation profile. If the 
confocal detection geometry is chosen properly a Gaus- 
sian effective volume with radius ro and height z$ is a 
good approximation |l8j |. 

In this work we want to focus on the possibility to 
measure flow velocities as it was shown by Gosch et al. 
0. 

The experimental autocorrelation function gi (t) is cal- 
culated from the measured fluorescence intensity F(t) by 



(F(t )F(t +t)) (5F(t )5F(t + t)) 
92(t) = (F(t )) 2 = 1 + <iW ' (4) 

For Brownian diffusing particles the concentration C(f, t) 
is known and can be used to calculate an expected auto- 
correlation function by combining equation and equa- 
tion © 



92(t) = l+-G DiS (t) 



(5) 



with N denoting the average number of marked objects 
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in the effective focai voiume and 

1 



TD 



t 



(6) 



where td = xjj denotes the mean passa ge time of 
the particle through the effective volume [2(J. and 
/ = zq/tq gives the aspect ratio of the effective volume. 
With the Rhodamine 6G diffusion constant Drqq = 
280 /im 2 s _1 [12j as a reference, the focal waist size and 
shape parameter were determined (ro — 0.20 (ira and 
/ ~ 8), and used for all subsequently acquired data. 

Equation © turns out to be the solution for the ideal 
case of point-like particles, which is well describing the 
case where the particle size R is small compared to the 
focal volume size r . For situations where R is compara- 
ble to, or even larger than ro no analytically closed form 
of equation © can be given (particle size effect). The 
autocorrelation function depends strongly on the distri- 
bution of fluorophores on the particle. For spherical par- 
ticles with radius R an approximate diffusion time can 
be given plj . 



R 2 



TD = 



AD 



(7) 



In the presence of a uniform translation in addition to 
the diffusive particle motion the autocorrelation function 
becomes 



where 



92 (t) = l + ^G Diff (i) GWw(f) 



G F iow(i) = exp( ) 



(8) 



(9) 



Taking into account the particle size effect the dwell time 
t f of the particles in the focus is given by 



ro + R 



(10) 



As can be seen from equation (JHJ the shape and the char- 
acteristic decay time of the autocorrelation function is 
dominated by the faster of the two competing processes. 
In order to obtain the flow velocities without pertur- 
bation by diffusion we used fluorescent microspheres of 
1 /i 77i radius (Fluosperes Yellow-Green 505/515, F-8823, 
Molecular Probes, Eugene, OR). On the experimentally 
relevant time-scale the Brownian motion of these spheres 
is of minor importance. Therefore we could neglect the 
diffusive term in equation ©. 

We used a commercial setup by Carl Zeiss (Jena, Ger- 
many) consisting of the module Confocor 2 and the mi- 
croscope model Axiovert 200 with a Zeiss C-Apochromat 
40x NA 1.2 water immersion objective. For excitation, 
the 514 nm line of a 25 mW Argon laser was used. Emit- 
ted fluorescence was detected at wavelengths longer than 
530 nm by means of an avalanche photo diode allowing 
for single-photon counting. 




FIG. 4: A dye jet pumped by the SAW transducer, shortly 
after switching the RF-power. In the middle of the transducer 
the fluid is pumped upwards, on both sides the fluid is sucked 
in. The bar is 0.5 mm wide. 




FIG. 5: A dye spot after 5 min of SAW pumping (40 mW) is 
shown. Superimposed is the edge line calculated by the image 
procedure. The bar is 1mm wide 



IV. RESULTS 

As a test for the image processing, a spot of Fluorescein 
pipetted into the hole of the cover plate was observed for 
10 min., with the SAW switched off. The area covered by 
the dye therefore only increases by diffusion. Using the 
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FIG. 6: The edges of a dye spot during the pumping of the 
SAW transducer (40 mW) are recorded every 10 s and put 
together in one graph. The evolution from a round spot to a 
butterfly shape can be observed. 



FIG. 7: The distance of the dye front as a function of time 
for 10, 25, 63 and 158 mW RF-power (from bottom to top). 
The data points where determined every 10 seconds at the 
line, marked by the arrows in Fig. [S] and connected by lines. 



linear diffusion equation (a; 2 ) = 2Dt, a value for the dif- 
fusion constant of the dye molecules of D=370/i m 2 / s can 
be determined which is in good agreement with literature 
values of 200-400>m 2 / s. 

In Fig. 0] a close up image of the area near the trans- 
ducer is shown, shortly after starting the SAW. The 1mm 
wide IDT is hidden by the dye solution on the lower side 
pumping a 200 /im wide jet to the upper side. The jet 
becomes broader with increasing distance and a part is 
transported to the influx areas of the IDT perpendicular 
to the direction of pumping. After 10 min., the Fluo- 
rescein covered area shows a butterfly shape. Employing 
the image processing the edge of the dye spot is deter- 
mined every 10 s during the observation time. A border 
line between the colored and uncolored areas computed 
this way is plotted in Fig. [5] Collecting the traces of a 
10 min. experiment into one graph yields a plot as in 
Fig. |SJ The initial dye spot and the succeeding evolu- 
tion of the edge lines are visible in this graph. This data 
was used for determining the evolution speed of the area 
(Fig. |8J . The growth of a Fluorescein spot was observed 
for 7 different RF-powers and the actual position of the 
front line in the direction of the arrows in Fig. was 
recorded. The growth speed decreases fast with distance 
from the transducer, until the diffusion constant of the 
dye (370 /J.m 2 / s) is reached. In the experiment, this po- 
sition x* D is found to depend linearly on the logarithm of 
the RF-power, see Fig.|5] 

From the movement of the dye front in the symmetry 
direction marked by the arrows in Fig.[S]the position can 
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FIG. 8: The speed of the growth of the area of the Fluorescein 
spot is plotted against the particular position of the advancing 
dye front for 4 different RF-power values (10, 40, 63 and 158 
mW from left to right). When the speed reaches the region 
of the diffusion constant of the dye, the data is superposed 
by noise. The diffusion constant of Fluorescein (370 fim 2 / s) 
is marked by a dashed line. For small distances, representing 
early times, the image processing yields misleading results, 
since the transducer region is covered by the initial round dye 
spot. 
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FIG. 9: The distance x* D where the speed of the area growth 
reaches the diffusion constant of Fluorescein (370 /xm 2 / s) de- 
pendening on the applied RF-power. Superimposed to the 
data is a straight line. 
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FIG. 10: The speed of the propagating dye front in the direc- 
tion of the arrows in Fig. |S| is plotted against the particular 
position for 4 different RF-power values (10, 40, 63 and 158 
mW from left to right). The speed values are calculated by 
differentiating the front position vs time data. Small dis- 
tances are not displayed, because the image processing yields 
misleading results close to the transducer. 
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FIG. 11: The speed of the advancing dye front at a distance 
4mm from the SAW transducer is plotted against the RF- 
power. A straight line has been fitted to the data. 



be determined as a function of time (Fig. 01 . Differentiat- 
ing this data and plotting it against the x-position of the 
front leads to a fast decrease of the speed (Fig. lTUj) . sim- 
ilar as found for the growth of the dye covered area. The 
decay of the velocity with the distance does not change 
qualitatively upon varying the RF-power. As a quantita- 
tive measure we show the velocity at the distance 4mm 
from the transducer in dependence of the applied RF- 
power in Fig. 1111 . In the measured range a linear depen- 
dence yields a reasonable fit although the extrapolation 
to small powers becomes unphysical. 

Fig.ll2lshows FCS correlation data taken as a function 
of the applied power. The correlation coefficient Gfiow(^) 
was determined from the measured intensity autocorrela- 
tion by means of the eqs. J3J and 10 under the assump- 
tion that Goiff(^) = 1) and compared to the exponential 
law @. The exponential fits well agree with the exper- 
imental data. The small deviations may be attributed 
to particle size effects [2lJ or to the polydispersity of the 
used beads [2i|. The fitted dwell times of the particles 
show a strong dependence on the applied power. Accord- 
ing to equation (|10f> the velocity was calculated and plot- 
ted as a function of the power (Fig. 1131) . This was done 
for three different positions in the channel. The highest 
velocities were achieved at positions near the chip. The 
measured data can reasonably well be fitted by a linear 
dependence of the velocity on the applied power. 

In Fig. E| the velocity profile for a selected area close 
to the chip is shown. The ^-direction is the direction in 
which the chip is supposed to pump. To investigate the 
dependence quantitatively we took the values indicated 
by the open squares and plotted their logarithm as a 
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FIG. 14: Velocity profile for a selected area near the pumping 
chip at a RF-power of 79 mW. The triangle, square and circle 
mark the positions at which the velocities shown in Fig. 1131 
were taken. 



FIG. 12: Normalized autocorrelation curves from FCS mea- 
surements. The applied power was varied between 150 mW 
(right) and 500 mW (left) in steps of 50 mW. Circles represent 
the experimental data whereas the lines are fits to equation 
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FIG. 13: Flow velocity as a function of the applied power 
for three different positions indicated by circles, triangles and 
squares. The corresponding positions are marked by the same 
symbols in the Fig |14l Circles represent the results from the 
fitted curves in Fig. IH The line is a fit to the data with a 
slope of 1. Flow velocities close to 10 fim/s are not determined 
correctly (points with large error bars). 




FIG. 15: Velocity along the central x-direction at the po- 
sitions shown by squares in Fig. 1141 Ideally, the velocity at 
the symmetry point x — should vanish. The high gradi- 
ent of the velocity field there and the experimental positional 
uncertainty, however, lead to the apparent finite velocity. 



function of the position, see Fig. ED As can be seen 
from this graph, the velocity decreases in a similar way 
as found for the propagation of the dye front in Fig. 
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THEORY 



Calculation of the flow field 



All observed flows in the 200/mi thick fluid film had 
velocities smaller than 1 mm/s and, hence, are charac- 
terized by Reynolds numbers smaller than 1 . This means 
that the flow is completely laminar |23j |. which simplifies 
the situation for the calculation of the streamlines. Al- 
though the detailed mechanism of the sound conversion 
into steady streaming is fairly complicated the steady 
flow is determined by the relatively simple Stokes equa- 
tions H : 



/jAv(x) = Vp(x) - P(x) 



V 







(11) 



(12) 



where fj, denotes the shear viscosity of the fluid. The 
driving force F is given in terms of the sound wave by 
eq. (J2J. The SAW force field is generated in the area of 
the transducer and decays exponentially with distance. 
At about 1 mm it has dropped more than three orders 
of magnitude. We use these features in order to find 
an approximate analytical solution of eqs. (flip and Q12|l 
for the time averaged velocity and pressure fields and 
simplify the problem by the following assumptions: (i) 
lateral boundaries are neglected, i.e. we consider an in- 
finitely extended fluid layer of thickness d with no-slip 
boundary conditions at both confining plates; (ii) we ap- 
proximate the actual spatially extended force density by 
a pair of constant forces being concentrated on lines per- 
pendicular to the plates separated by the distance a: 

F(x) = (f+5(x - a/2) + f_5(x + a/2)) S(y) (13) 

where (x, y) = r are coordinates spanning the planes par- 
allel to the plates and z is a coordinate perpendicular to 
the plates with z — on lower and z = d on the upper 
plate. The two forces f± are assumed to have the same 
component in the z-direction and opposite components 
in the direction of the vector separating the two forces, 
i.e. in the x-direction, pushing the fluid away, 



f± = f z e z ± fe x 



(14) 



where e XiZ are the unit vectors in x- and z-direction, re- 
spectively. We expect that outside the narrow spatial do- 
main where the actual force does not vanish this presents 
a fairly good approximation to the actual situation. Be- 
cause of the linearity of the Stokes equations (|11I12|> it 
is sufficient to consider a single line-force contribution. 
The final result can be obtained by the superposition of 
the solutions to the two forces. 

Having only one line-force we may put the position of 
coordinate center at the position of the force. Further 
we note that one formally may include the z-component 



of the force into the pressure term by introducing the 
effective pressure 



p eff (x) = p 2 (x) - f z zS(x)S(y) 



(15) 



Measuring lengths, velocities and pressure in units of d, 
f/fi and f/d, respectively, leads to the dimensionless 
Stokes equations: 



Av(x) - Vp cff (x) 
V • v(x) = 



-e x S 2 (r) 



with no-slip boundary conditions on the plates: 



v(r,0) =v(r,l) 







(16) 
(17) 



(18) 



where S 2 (r) = S(x)S(y) is a two-dimensional <5-function. 
If one takes the divergence of the Stokes equation one 
finds a Poisson equation for the effective pressure: 



Ap cff (x) = |-£ 2 (r) 



(19) 



Because the inhomogeneity on the right hand side is z- 
independent it follows that the effective pressure also 
does not depend on z. Consequently the z-component 
of the velocity fulfills the Laplace equation with homoge- 
neous boundary conditions and therefore v z {x) vanishes 
everywhere. The remaining two components of the veloc- 
ity field can be expressed in terms of a stream function 
$(x): 



Ux(x) 



<9$(x) 
dy 



<9$(x) 

dx 



(20) 



Differentiating the x- and y-components of the Stokes 
equation with respect to y and x, respectively, and sub- 
tracting them one obtains for the stream function: 



A r A$(x) = -|- ( 5 2 (r) 
ay 



(21) 



where A r = d 2 /dx 2 + d 2 /dy 2 denotes the two- 
dimensional Laplace operator in the (x, y)-plane. Bound- 
ary conditions for $(x) are: 



$(r,0) = $(r,l) =0 
<9$(x) _ <9$(x 



dx dy 
lim $(r, z) = 



= 



(22) 
(23) 
(24) 



where r = (x 2 + y 2 ) 1 ^ 2 denotes the absolute value of r. 
Because of the periodic boundary conditions with respect 
to z the stream function can be represented as a sin- 
series: 



$(x) = c„(r) sin(7rnz) 



(25) 
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Here the coefficients are still functions of r and obey the 
equations 



A r (A, 



(im) 



where 



In 



2 / dz sin(7rn.z) 
lo 





4 

nil 



for n even 
for n odd. 



(26) 



(27) 



The solutions of these equations can be expressed in 
terms of the Green's functions of the time-independent 
two-dimensional Klein-Gordon equation: 



(A r -£: 2 ) ffm (r) = 5 2 {v) 



(28) 



where k is a positive number. This Green's function is 
given by |24|: 



9k(r) = -—K (kr) 
Zir 



(29) 



where Ky(z) is the modified Bessel function of order u, 
see Ref. [25j • To achieve this goal we introduce an auxil- 
iary function and rewrite eq. 126[l in terms of two Klein- 
Gordon equations: 

A r Mr) = -|"<5 2 (r) (30) 
ay 

(A r - (n7r) 2 )c„(r) = 7 „/i(r) (31) 

The solution to the first equation can readily be expressed 
in terms of the Green's function for m = 0: 

d 

h(r) = - lim ^-g m (r) 

m— >oo ay 

i a 

= — — In r 

2tt ay 

From the second equation one then finds: 
c„(r) = cZ 2 r'. 9mr (r-r')M r ') 



(32) 



In 9 

(2tt) 2 9y 



J (fr' K (nnr')ln\r-r'\ (33) 



Introducing polar coordinates one may perform the inte- 
gration over r and obtains: 



c„(r) = 







for n even 



^ 2 (1 — niTrKi(rnTr)) for n odd 



(34) 



The modified Bessel function decays rapidly for large val- 
ues of its argument and therefore nnrKi (mrr) can be ne- 
glected for 7rr 3> 1. This gives the asymptotic behavior 
of the stream function: 



$ as (x) 



V 



z(l-z) 



(35) 



where we summed the smc-senes > „ „j j — v \ ' = 

L — *n oaa n 

^j-z(l — z), see Ref. [2(|. Deviations from this asymp- 
totic result can only be seen within a distance of a few 



$(0,y, 0.5) Q 



02 



0.01 - 



| i i i i | i i i i | 




i i i i | i i i i | i i i i | 
1 2 3 

- -0.01 



0.02 



FIG. 16: The streaming function $(0, y, 0.5) (solid line) per- 
pendicular to the axis x — displays extrema close to 
\y\ — 0.3. The broken line shows the asymptotic far field 
behavior according to eq. 13511 . Both functions perfectly co- 
incide for distances from the origin that are larger than twice 
the plate distance. 




FIG. 17: The streaming function $(0, y, z)/$(0, y, 0.5) di- 
vided by its maximal value at z = 0.5 at x = is shown 
as a function of the transversal coordinate z for different dis- 
tances y from the origin where the force is sitting. In the 
very close vicinity of the force source, this profile is flat in the 
middle part of the fluid layer. For larger distances from the 
force it rapidly approaches the asymptotic parabolic profile, 
see eq. I|35|l . 



gap widths from the location where the force is acting, 
see Figs. Ijlfijl and i|17)l. The streaming function of the 
force dipole readily follows as the linear combination of 
two velocity fields that are caused by slightly displaced 
line-forces pointing in opposite lateral directions: 

$di P oic(x) = *(x - ee/2) - $(x + ee/2) (36) 
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t = 15 



t = 50 



t = 100 
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FIG. 18: The graph shows the calculated evolution of a dye 
spot pumped by a SAW transducer in the center of the images. 
The black arrows mark the direction of the velocity field. The 
initial dye spot is marked by the circle. The 4 black ellipses 
in the lower image represent areas where the dye is trapped 
by an eddy. The trajectories were calculated for the plane in 
the middle of the fluid film with diffusion neglected. 



where e is the width and e the unit vector in the direc- 
tion of the displacement. In the following we will assume 
that both the displacement and the forces show in the 
x-direction. Using the asymptotic expression (|35[l and 
taking the leading contribution with respect to the dis- 
tance e one obtains: 



dipole 8?rr 4 



z(l-z) 



(37) 



The velocity field then becomes asymptotically: 



B. Motion of advected particles 

We consider the motion of a small particle that is ad- 
vected by the flow field iJSHJl. Its trajectory agrees with 
the streamlines of the vector field v(x) which are the 
solutions of 



x(t) = v(x(t)) 



(39) 



Because the velocity component in the z direction van- 
ishes the motion of a particle is strictly confined to the 
initial plane z = z(0) as long as diffusion can be ne- 
glected. Fig. 1181 shows the form that an initially circular 
spot in the central plane (z = 1/2) acquires at four later 
times. We note that the area covered by the spot does 
not change because the vector field has vanishing diver- 
gence. The apparent change of the area of a dye spot 
which is observed in the experiment is due to the fact 
that different layers that move differently are superim- 
posed. The alternating stripes of colored and uncolored 
fluid that evolve after sufficiently long time within each 
vertical layer cover each other and therefore cannot eas- 
ily be seen in the experiment. In the course of time these 
layers become ever thinner and eventually diffusion be- 
comes effective. Around the center of each of the four 
eddies regions exist that remain invariant under the flow, 
i.e. there is no transport out of these regions or into them 
other than diffusion. 

Finally we consider a particle that initially sits in the 
middle position between the plates on the x-axis at a dis- 
tance xq from the origin. It moves with the fluid strictly 
in the x-direction according to 



dt 



x(t) = v x {x(t), 0,0.5) 



x(tf 



(40) 



where c = e/32w is a constant. The solution of this dif- 
ferential equation reads: 



x{t) 



ct 



1/4 



(41) 



™ (i_ 4 »j ]*(!-*) 
iirr \ r z 



QJ 



87rr 4 



I-4-3 2 (1- z) 



(38) 



Note that the near field has a less pronounced, 
logarithmic singularity. The full two-dimensional 
vector-field v x (r, 0.5), v y (r, 0.5)) that follows from the 
cas. H2QI25l34l3(j|) is shown for a pair of opposite line 
force pointing in the positive and negative ^-direction 
at r = (±1,0) in Fig^| In this case the fluid is sucked 
along the y-axis and pushed out in the ±ir-direction gen- 
erating an eddy in each quadrant close to the dipole. 



Here we have neglected the influence of diffusion. This 
will be of primary importance for the motion in the z 
direction because there is no systematic force in this di- 
rection and, moreover, the z-coordinate of the particle 
position strongly influences its lateral velocity. A dif- 
fusional motion in the z-direction averages over z and 
results in an effective value of the constant c which be- 
comes smaller by the factor 2/3. 

In order to compare with the theoretical prediction of 
eq. (|41|l in Fig. ^3 the fourth power of the position is 
plotted against time for four different RF-power values. 
Except for the highest power value where the asymptotic 
velocity field eq. (|4Tj)l might not be fully applicable, the 
experimental data nicely fall onto straight lines. From 
the slope of this line we obtain the constant c which 
determines the velocity field, eq. (|40|l . The value of c 
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FIG. 19: The fourth power of the front positions (solid line) 
for four values of RF-power (10, 40, 63, 158 mW from bottom 
to top) of the dye spot indicated by the arrow in Fig. |S| mea- 
sured at different instants of time (see also from Fig. |7J falls 
on straight lines as predicted by the asymptotic law given 
in eq. 14111 . Only for the very first points and at high RF- 
power deviations appear. From the inclination of the linear 
fits (broken lines) one obtains a reliable fit of the constant c 
that determines the velocity field, see eq. 14011 . 

increases with the applied RF-power, see Fig.l2UI An ex- 
trapolation of a linear fit of the data, however, would lead 
to a spurious finite positive value of c and consequently 
to finite velocities in the absence of driving. We therefore 
conclude that the increase of c with the RF-power must 
be weaker than linear. 



VI. DISCUSSION 

We applied fluorescence microscopy and FCS for the 
flow profile analysis of a SAW transducer in a 200/im 
thin fluid film. With the first technique the growth of 
the area of a Fluorescein spot was observed. From this 
the flow speed in x-direction could be determined. 

In the FCS experiments, the flow speed in a homoge- 
neous dispersion of fluorescent beads was measured at 
different points in a 5mm wide channel. With this tech- 
nique no limitations exist for the choice of the position 
where the data is collected. 

Both methods yield a fast decay of the flow speed with 
distance from the transducer with a comparable rate. As 
a theoretical model, the Navier Stokes equations were 
solved for an infinite fluid layer in the Stokes limit of low 
Reynolds numbers. The conversion of the sound wave 
into a steady flow was described in terms of body forces 
that were further simplified as line forces. Because the 
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FIG. 20: The constant c as determined from measurements 
of the position of the dye front is shown for different values 
of the RF-power. 



actual spatial extension of the body forces is confined to 
a narrow regime this seems to be a reasonable approx- 
imation. The resulting flow field is in good qualitative 
agreement with the experimental findings. According to 
the theoretical prediction, the velocity field asymptoti- 
cally decays with the third inverse power of the distance 
from the IDT. The resulting motion of the front of a dye 
spot in the symmetry direction is in quantitative agree- 
ment with the measurements. 

To determine the distance at which the influence of 
the transducer on the dye molecules vanishes, the growth 
speed (va) of the fluorescein covered area was measured. 
Plotting this data against the actual position of the ad- 
vancing front shows also a fast decay until the diffusion 
constant of the the dye is reached. This distance depends 
logarithmically on the RF-power and is about 8 mm for 
the maximum applied 158 mW in the microscope study. 

When the flow speed is measured at different points 
in the streamlines with both methods, a linear behavior 
dependent on the RF-power can be observed. However, 
when the constant c in eq. I|41|l is determined from the 
x position vs time data for different Hf-values, no line 
through the origin can be fitted. So following eq. (|40|l the 
flow speed cannot depend linearly on the power at a cer- 
tain point in the flow field. The power range where mea- 
surements where taken, therefore probably represents an 
approximately linear regime. 

Finally we address the question what one can learn 
from the presented experiments and theory in order to 
accelerate the hybridization of DNA in a thin fluid film 
with the lateral size of a normal glass slide (75x25mm). 
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For this size, diffusive transport becomes extremely slow 
and should be complemented by advected transport. The 
flow induced by surface acoustic waves can be efficiently 
used for this purpose. However, because of the rapid 
decay of the induced velocity field away from the trans- 
ducer the maximal mixing range was found to have a 
linear dimension of about 8 mm in the present experi- 
ment at the highest RF-power. A further increase of the 
power does not seem feasible because it would enlarge 
the mixing range only insignificantly. To account for the 
limited range of the induced flow by the SAW transducer 
the number of transducers on a slide should be increased. 
In the calculated flow field in Fig. ^] four eddies appear 
close to the transducer where particles can be trapped. 



To avoid this problem, pairs of transducers should al- 
ways be operated alternately. So the trapped particles 
from one transducer will be moved by the other and vice 
versa. 
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